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We construct a multimetric gravity theory containing N > 3 copies of standard model 
matter and a corresponding number of metrics. In the Newtonian limit, this theory generates 
attractive gravitational forces within each matter sector, and repulsive forces of the same 
strength between matter from different sectors. This result demonstrates that the recently 
proven no-go theorem that forbids gravity theories of this type in N — 2 cannot be extended 
beyond the bimetric case. We apply our theory to cosmology and show that the repulsion 
between different types of matter may induce the observed accelerating expansion of the 
universe. In this way dark energy can be explained simply by dark copies of the well- 



o 

CN understood standard model. 



I. MOTIVATION 



The widely accepted standard model in modern cosmology is known as the ACDM model. 
Its theoretical basis is a homogeneous and isotropic spacetime metric with dynamics governed by 
general relativity. Already this simple setting allows for a successful explanation of very different 
astronomical observations, such as the cosmic microwave background [I], the accelerating expansion 
of the universe [21 [3] , and its large scale structure [I] . This explanation requires that the visible 
£f~) standard model matter only contributes about 5% to the total matter content of the universe 

o 

and must be augmented by an incredible 95% of dark matter and dark energy. However, the 
constituents of dark matter and dark energy are not specified by the ACDM model, and their 
^ nature is presently unknown. 

This situation has led to the development of numerous models for dark matter and dark energy, 
both from the perspectives of particle physics and of gravity. Particle physics models for dark 
matter [5] include weakly interacting massive particles [6], axions [7], or massive compact halo 
objects [SJ. Dark energy [9j is modelled e.g. by scalar fields as quintessence [TUl HTj or K-essence P21 
H5] , as a Chaplygin gas [Hj, or by employing tachyons. In contrast to these particle theoretic 
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approaches, modifications of general relativity may be employed in order to explain the effects 
which are otherwise attributed to dark matter or dark energy. The simplest example of such 
a modification is the introduction of a cosmological constant. Other examples include modified 
Newtonian dynamics |15j , tensor vector scalar theories \1Q\ [T7] , curvature corrections by the full 
Riemann tensor as in [18] or by the Ricci scalar in f(R) theories [19], higher-dimensional models 
such as the DGP model |20l [2"T] , or structural extensions such as non-symmetric gravity theory [52] 
and area metric gravity |23|, [2^] . 

In a recent article [25] we speculated on the possibility that both dark matter and dark energy 
might be constituted by an additional copy of the standard model which couples to a second metric 
and interacts with visible matter exclusively through gravitation. If the gravitational coupling is 
attractive within each matter sector, but repulsive between the two different matter sectors, this 
might explain the dark universe. Moreover, matter from a dark copy of the standard model could 
be located in the so-called galactic voids, which are seemingly empty if one considers only visible 
matter. Their presence then would assert a repulsive force acting on visible galaxies, so that these 
are pushed away from the galactic voids. Astronomical observations indeed suggest the existence 
of forces of this type |26| I27j . The major advantage of such a model would be its clear and 
simple interpretation: since dark matter and dark energy are simply constituted by a dark copy 
of standard model matter, their physical properties would be well-understood. This model would 
obey a version of the Copernican principle: it seems more natural to assume that the dark universe 
is constituted by the same type of matter known from the visible universe, than that the visible 
universe is distinguished from the dark universe by its physical properties. 

Unfortunately, gravity theories of the type on which this speculation is based are not so easily 
realized. We could prove a no-go theorem in [25] that exludes all canonical bimetric extensions 
of Einstein gravity, i.e., all those bimetric theories with two copies of standard model matter that 
have a Newtonian limit in which the attractive gravitational forces within each matter sector and 
the repulsive forces between matter belonging to different sectors are of equal strength. 

In this article we will demonstrate the existence of canonical extensions of Einstein gravity with 
attractive and repulsive forces beyond the bimetric case, i.e., for N > 3 metrics and a corresponding 
number of standard model copies. In section [IlJ we will explicitly construct an action for such a 
theory, derive its equations of motion, and show that the theory has the required properties. We 



will then discuss the cosmological consequences of our theory in section III For a simple model we 
will show that the accelerating expansion of the universe conventionally attributed to dark energy 
can indeed be induced by the mutual repulsion between several copies of standard model matter. 
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We will conclude with a discussion in section HVl 

II. CANONICAL EXTENSION OF EINSTEIN GRAVITY 

In this section we will show the existence of gravity theories with N > 3 metrics and a corre- 
sponding number of standard model copies which have a Newtonian limit in which the attractive 
gravitational forces within each matter sector and the repulsive forces between matter belonging 
to different sectors are of equal strength. Starting from a set of well-motivated assumptions, for in- 
stance that all but one sector appear dark for any observer, we will first construct an action ansatz 
for such a theory. Second, we will derive the corresponding equations of motion by variation and 
show that the remaining parameters in the action can be chosen so that the theory acquires the 
proposed Newtonian limit. In the following section we will analyze the cosmological consequences 
of the theory. 

A. Construction of the action 

The basis for our extension of Einstein gravity is a four-dimensional manifold. The field content 
we consider is given by a set of N > 3 metric tensors g 1 , . . . , g N and also by N copies of the standard 
model with fields which are collectively denoted by . . . , ^> N . This reflects our motivation that 
dark matter and dark energy should be constituted purely by additional copies of the standard 
model without introducing any other field. A standard model copy together with the corresponding 
metric g 1 will be called a sector of our theory. To proceed towards an ansatz for a suitable action 
we will use the following assumptions: 

(i) The fields SB 1 of each copy of the standard model couple only to the corresponding metric g 1 . 

This assumption is needed in order to obtain the correct behaviour of matter within a gravi- 
tational field. The fact that each type of matter is affected only by a single metric guarantees 
that the motion of observers in the sector with metric g 1 is governed by the corresponding 
set of timelike geodesies, and the standard notion of causality of matter fields \P J is provided 
by the Lorentzian cones. 

(ii) Different sectors couple only through the gravitational interaction between the metrics. 

Since there is no non-gravitational evidence for the existence of additional standard model 
copies, we must assume that there is no direct non-gravitational coupling between them. In 
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consequence, matter from any given sector will appear dark for observers in all other sectors. 

(Hi) The equations of motion contain at most second derivatives of the metrics. 

This assumption is one of mathematical simplicity, and guarantees a reasonable amount of 
technical control over the partial differential field equations. It will be useful to restrict the 
possible terms in the action of our theory. 

(iv) The theory is symmetric with respect to an arbitrary permutation of the sectors (g 1 ,^ 1 ). 

This assumption is made for simplicity; it employs the Copernican principle in the sense that 
the same laws of nature should hold within each sector. It also follows that the interaction 
between the different sectors will satisfy Newton's principle that action equals reaction for 
the gravitational forces. 

Establishing assumption means that the action we look for must contain in its matter part 
a sum over copies of standard model actions 

where to 1 = d^x^/g^ denotes the canonical volume form related to g 1 , and CmW is the stan- 
dard model scalar Lagrangian. Assumption then implies that the remaining gravitational part 
of the action can only depend on the different metrics. Hence the total action can be decomposed 
in the form 

N 

S = S G [g\...,g N } + Y,SM[g 1 ,* 1 }. (2) 

7=1 

We now turn our focus to the gravitational part of this action, which can be written as 

Sg[ 9 \..-,9 N ] = 1 j mC G [g\...,g N } (3) 

for a symmetric volume form luq = d^x^/go with 

N 

90 = II (9 1 ) N (4) 
i=i 

and a scalar Lagrangian Coig 1 , ■ ■ ■ ,g N ]- We use units so that the Newton constant is normalized 



as 8ttGi\[ = 1 and [Co] = L~ 2 . As a consequence of assumption (Hi), the Lagrangian cannot 
contain terms with higher than second derivatives of any metric g 1 . Hence the only tensors that 
may appear in this Lagrangian are the metrics g 1 , the connection difference tensors 
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their covariant derivatives S7pS JK l jk, and the Riemann curvature tensors R Il jki for each metric. 
For simplicity, and in analogy to Einstein gravity, we now construct our Lagrangian only from 
terms of the form g Il ^Rj-. From assumption ^kfy we then deduce that the prefactor of each of 
these terms should be independent of the individual sectors /, J. However, it may still depend on 
whether I and J are equal or not. We therefore choose the following ansatz for the Lagrangian, 

N 

C G [g\ . . . ,g N ] = (* + yS IJ )9 Hj Ri ■ (6) 
Z,J=1 

The parameters x, y are constant, and imply that g 1% ^Rf^ appears with prefactor x + y if I = J 
and prefactor x otherwise. 

Equations Q and ^ define the gravity theory we wish to investigate in the following. 

B. Derivation of the field equations 

We will now derive the gravitational field equations from our action ansatz @ and © by 
variation. In particular we will show that the parameters x, y can be determined so that the theory 
obtains a Newtonian limit in which the attractive gravitational forces within each matter sector and 
the repulsive forces exerted from matter belonging to different dark sectors are of equal strength. 

The variation of the gravitational part of the action can be written in the form 

SSa = \ J2( x + ^ /J ) ( / d 4 ^vW^'4 + / M^Rfj + g^SRfj)) . (7) 
i,j=i \J J J 

It is straightforward to compute the variations of the occurring terms. We will therefore only give 
a brief sketch of the computation. For the variation of the volume form, note that 

N 



Vso = ^E^^. (8) 
1=1 

The variation of the inverse metrics is given by the standard formula 5g Il i = —g Ita g I: ' b Sg^ b . For 
the variation of the Ricci tensors, we use the formula 



5R{ 3 = (g^Sp^ - lg Ja %^ - lg Jc %S b ^ V J d V J c 5g J ab . 



(9) 



The occurring covariant derivatives on Sg^ can be resolved by repeated use of the partial integration 
formula 

"w VfV < = - / uo&iV*, (10) 
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which holds for arbitrary vector fields V. Here and in the following calculation we use a convenient 
short notation for contracted connection differences, 



s IJ %, (11) 



and for the arithmetic mean with respect to the first sector index 

N , N 



b i k = N ^ jkl * = N 2-" i ' ^ ' 

i=i i=i 



Further, note that covariant derivatives on the metrics can be written as V^g^, = —2S IJd a ^g^ id , 
using the fact that g J is covariantly constant with respect to V J . Thus we finally obtain the 
variation of the gravitational part of the action in the form 

N 

SS G = -^J2 ^ Iab Sgi b (13) 

7 1 J 



1=1 

with 

N X 

n Iab 

2N 



K Iab = ~^9 Iab £ (s + ZV^ + Hj + y^/V 

J,K=1 J=l 

- 2 (g'^Sp^ - lg Ia %S^ - \g Icd S^ J> + y 5 IJ ) fa^S^^S 1 d) (14) 

+ -g^&J^ + ^Vi&i + VlS IJi dp g J i p + S IJp cq S IM dp g J ™ + S'^S"' dp g Jpq ) . 

We still have to compute the variation of the matter part of the action with respect to the 
metric tensors. Since each type of matter couples only to a single metric tensor, this variation can 
be written in standard fashion in terms of the matter energy momentum tensors, 

5S M W,¥]= l -J u/t'^V (15) 

Note that this integral is performed using the volume form u 1 , whereas the variation SSq above of 
the gravitational part of the action contains the symmetric volume form ujq. This can be accounted 
for easily by recalling that ujq = oj 1 y^go/g 1 ■ Thus, by combining (13) and (15) we finally obtain 
the full equations of motion 



T Iab = ^g /g I K lab . (16) 

Note that the only maximally symmetric vacuum solution of these equations is g 1 = \ Ir q for 
constants A and flat Lorentzian metric rj. So the Newtonian limit of the equations can be obtained 
by linear gauge-invariant perturbation theory with the ansatz g 1 = X 1 (ry + h 1 ) where one assumes 
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small components |/i 7 b | -C 1. The computation can be performed in complete analogy to the 
bimetric case which is shown in detail in [25]. One needs to determine the dependence of the 
gauge-invariant Newtonian potentials & on the matter densities p 1 (whose definition absorbs the 
constants A 7 ). This calculation then results in Poisson equations 

N 

2 



A$ 7 



^E^V (17) 



j=i 



with a constant coupling matrix A IJ . Linearizing the equations (16) of our theory here, we obtain 
the matrix components 

3 V yj \4N(Nx + y) J v ' 

As discussed at the beginning of this section, canonical extensions of Einstein gravity in our sense 
are defined by a standard Newtonian limit within each matter sector; this is achieved by diagonal 
entries A n = 1 since 8ttGn = 1- Moreover, the canonical extensions have repulsive gravitational 
forces of equal strength between matter from different sectors, i.e., off-diagonal entries A IJ = — 1 
for I 7^ J. These two requirements are met for parameter values 

2N-1 -2N + 7 

X 6N(2 — N) ' V 6(2 - TV) ' 1 ' 

There are two immediate special cases. For N = 1 the action of our theory, see Q and (|6]), reduces 
to the Einstein-Hilbert action, and we obtain the Einstein equations with only one matter sector 
and standard Newtonian limit. For N = 2 the parameters x,y above are not defined, and this is 
consistent with our no-go theorem [25] for bimetric gravities of this type. Finally, for N > 3, this 
result verifies our proposition on the existence of canonical extensions of Einstein gravity. 



III. COSMOLOGICAL MODEL WITH ACCELERATING EXPANSION 

The construction of the previous section provides us with an explicit gravity theory including 
dark sectors and repulsive forces. We will now analyze some of its cosmological consequences under 
the standard assumption of a homogeneous and isotropic universe. We will argue that the very 
early and the very late universe should be amenable to an effective metric description where the 
metrics from all sectors have an approximately identical evolution. For this case we will compute 
the reduced equations of motion. We will show that these resemble the Einstein equations, except 
for an additional negative factor that depends on the number N > 3 of sectors and rescales the 
gravitational constant. From the cosmological equations of motion we will then read off several 



s 



features of our model: the universe must be open and its expansion is accelerating. We confirm 
this also by obtaining all explicit solutions for radiation and dust matter. The early universe turns 
out to feature a big bounce rather than a big bang, while the acceleration of the late universe 
naturally becomes small. 



A. Effective metric ansatz 



The extrapolation of the Hubble expansion of the universe back in time suggests that the early 
universe becomes increasingly dense and hot. All matter hence moves relativistically so that one 
may describe this early stage dominantly by radiation. In our multimetric theory we simply extend 
this assumption to all matter sectors. This is another instance of the Copernican principle which 
suggests symmetry between the different sectors. On this philosophical basis, it seems reasonable 
to assume that the initial conditions for all matter sectors were the same at some early time. The 
field equations that are symmetric under permutation of the sectors then allow for a very similar 
non symmetry-breaking evolution of the sectors. This common evolution applies both to the energy 
momentum tensors as well as to the metrics, if averaged over cosmological scales. By this argument 
the very early universe can be described by means of a single effective metric g^ b = g a b and single 
effective energy momentum T^ b = T a b- This argument will break down as soon as perturbations 
in the sectors start to grow. These will lead to local symmetry breaking which should eventually 
transfer to a different cosmological evolution of the sectors. 

The symmetry between all sectors in our theory suggests that we should extrapolate our ob- 
servation of the Hubble expansion to all matter sectors. At very late times the universe hence 
will have sufficiently expanded so that the matter in all sectors can be described by dust. This 
implies that the structure formed at an intermediate age of the universe no longer influences the 
cosmological evolution. Since the physical laws in all sectors are the same and the initial conditions 
at some early time agree, as argued above, it is a plausible assumption that the intermediately 
different evolution of the sectors averages out again at very late times so that the effective metric 
and effective energy momentum solutions of our theory become attr actors. 

With these arguments the effective metric description, g^ b = g a ^ and T* b = T^, becomes 
available as a simple model both for the very early and the very late universe. We will now discuss 



the consequences of this assumption under which the equations of motion (16) greatly simplify. 
First, note that the connection difference tensors S IJ defined in ([5]) all vanish; this is due to the 
fact that all connections are equal. Second, the sums over sector indices /, J in the equations of 
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motion can be performed explicitly. This procedure results in precisely the same effective equation 
for each of the sectors: 

(2-N)T ab = R ab -^Rg ab . (20) 

Remarkably, these are the Einstein equations except for the additional factor 2 — N, which acts as a 
rescaling of the gravitational constant. Our assumption N > 3 implies that this factor is negative. 
Thus, the sign of the gravitational constant flips, and gravity for the effective metric becomes 
repulsive. By specializing the effective metric to the Robertson-Walker form in the following 
subsection, we will see that this results in an accelerating universe. 

B. Cosmological equations 

Homogeneous and isotropic cosmologies are characterized by the existence of six Killing vector 
fields responsible for spatial translations and rotations. The requirement that these fields are 
symmetry generators for all metrics g 1 in our theory restricts their form to be of Robertson- Walker 

type, 

g 1 = -nj(t)dt ®dt + ajit^^dx <g> dx 13 , (21) 

with lapse functions ni (t) , scale factors aj (t) , and a common purely spatial metric "y a /3 of constant 
curvature k £ {—1, 0, 1} and Riemann tensor R(l) a p-y5 = 2&7a[775]/3- Note that the lapse function 
ni in a single-metric theory may be set to unity by an appropriate rescaling of the cosmological 
time t. In a multimetric theory, however, there are N independent functions nj, which cannot be 
set to unity simultaneously. 

The matter content consistent with the cosmological symmetries is given by a set of N homoge- 
neous fluids with density pj(t) and pressure pi(t). Their energy- momentum tensors can be written 
as 

T Iab = ( P/ + Pl )u Ia u Ib + Pl g Iab (22) 

with velocities normalized by the relevant metrics from their sector so that g^ b u Ia u Ib = — 1. These 
tensors can be decomposed into the components Tq = pjnj and T^g = pia 2 j^f a p. 

We will now restrict this general multimetric cosmological model to the effective metric case 
of relevance to the very early and very late universe, as argued in the previous subsection. This 
means we can omit the sector index I from all matter densities and pressure functions as well 
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as from the lapse functions and scale factors. We may now rescale the cosmological time so that 



n(t) = 1. Since all metrics are now equal, we can use the symmetric field equations (20) and insert 
the Robertson- Walker metric ( |21[ ) (without sector index) . This leads to the cosmological equations 
of motion wherin dots denote derivatives with respect to t: 



3 (a 2 k , 



1 (a a 2 k . 



2 - N V a a 2 a 2 
The second equation can be replaced equivalently by the continuity equation 

p=-3^(p + p), (24) 

which can be derived alternatively from energy momentum conservation V a T a0 = which is a 
consequence of diffeomorphism invar iance. 

The first crucial observation from the above equations is that the matter density p can only 
be positive, if the universe is open with k = — 1. This is a prediction of our simple cosmological 
model, and contrasts general relativity where cosmological solutions for all three cases k = 1, 0, — 1 
exist. We further see that positive p constrains a by the inequality a 2 < 1. 

Without solving the equations of our cosmological model, we may obtain another remarkable 



result: an accelerating universe. To see this we form a suitable linear combination of equations (23) 
to obtain the acceleration equation 

a N -2 
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(p + 3p). (25) 



The strong energy condition, 



*- \Tg ab y a t b >Q (26) 

for all timelike vector fields t a , holds for all standard model matter. For perfect fluid energy 
momentum and using t a = u a this implies p+3p > 0. Since N > 3, it then immediately follows from 
the acceleration equation that a must be positive. This is a major difference to the cosmological 
solutions obtained in Einstein gravity with N = 1, where a positive acceleration cannot be obtained 
without either a cosmological constant, or an exotic type of matter which has sufficiently negative 
pressure p < — p/3. Within our theory, the acceleration is caused solely by the fact that the sign 
of the gravitational constant which is sign(iV — 2) flips for N > 3. 
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C. Explicit solution 

We will now find the exact solutions to the cosmological equations of motion of our model. 
These will explicitly confirm that the accelerating universe also is expanding. We will see that the 
acceleration tends to zero for very late times, and that the early universe features a big bounce, 
not a big bang as it does in Einstein gravity. 



In order to solve the cosmological equations (23), we introduce the conformal time parameter rj 



which is defined by dt = adr). Denoting derivatives with respect to rj by a prime ', we obtain the 
open universe k = —1 equations 

N V a 4 a 2 
1 f„a" a' 2 1 



(27a) 



2 — TV \ a a a 4 ar / 
Applying these equations to the radiation-filled early universe requires an equation of state p = ojp 
with equation of state parameter oj = 1/3, while the late universe requires the choice u = for 
dust matter. Inserting the equations of motion above into the equation of state leads to 

= U p - p = (2 _ 1 jV)fl4 {2a" a + (3w - l)a' 2 - (3w + l)a 2 ) . (28) 

The general solution of this equation takes the form 

/3w + l \ / 3w + l 
oi ex P I — 2 — V j + a 2 ex P ( 2 — V j j ' ' 

for integration constants a\,a%. Employing this explicit expression for the scale factor in equa- 



tion (27a) we compute the matter density 



12 / ^ + 1 \ / 3^ + 1 \Y%£ 
P= N _ 2 aiCl2 ( Q i ex P( — 2 — ^)+ a 2 ex P( 2 — V j) ' ^ ' 

The values of the integration constants in this solution are constrained by the requirement 
that a and p should be positive. This can only be achieved if both a\ and 02 are positive. Then 
it is not difficult to check another important feature of the solution: the scale factor a attains a 
positive minimal value 

a = (Aa^)^ 1 (31) 

at conformal time 

Vo = In ^ . (32) 

3co + 1 a% 
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This property of our cosmological model tells us that every solution features a big bounce where 
the matter density becomes maximal, 



3 

(N~2)al 



PO = 777 7772 • ( 33 ) 



An alternative way to parametrize the solutions uses the values rjo and ao at the big bounce 
instead of the original integration constants a\ and a2, which yields 



a = ao 



( /3w + l \\37+i 

(coshf — - — (v-Vo)JJ , (34a) 

6^+6 

/ /3o; + l \\ 3^+1 

p = poi coshf — - — (v-Vo)JJ ■ (34b) 

From this representation of the solutions one immediately sees why ao and po are extrema of the 



scale factor and matter density, respectively. Using the definition dt = adrj and (34a) we may 
transform the solutions back to cosmological time. For general oj the integrated relation between 
t and 77 is 

t = 2*1 ( -^-i -e^+^-ri) , (35) 

43^+1 e v-vo ^3w + l'3w + l 3w + l / 

in terms of the hypergeometric function 2-F1. Now the big bounce at 77 = rjn corresponds to t = 0. 




FIG. 1: This figure shows the scale factors of the radiation-filled universe (dashed line) and the dust- filled 
universe (solid line) plotted over cosmological time. 

The early universe (near the big bounce) is filled with radiation with u = 1/3. In this case the 
solutions simply become 

± = y/l + t*/4, ^ = {l + f/alY 2 . (36) 

an v On 
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These are plotted as the dashed lines in firgures [T] and [2] The late universe is modelled by dust 
with oj = 0. For this case we may consider the asymptotic behaviour of the acceleration. One can 
check that the acceleration a as a function of rj is 

(37) 



a (1 + cosh(ry - 7] )) 2 

For late times t — > oo which correspond to ry — > oo the acceleration tends to zero. The dust 
solutions for the scale factor and matter density are plotted as solid lines in figures [T] and [2] 




FIG. 2: This figure shows the matter densities of the radiation-filled universe (dashed line) and the dust-filled 
universe (solid line) plotted over cosmological time. 



IV. CONCLUSION 

In this article, we have presented a multimetric extension of Einstein gravity with the field 
content of N > 3 metric tensors and a corresponding number of standard model copies. The 
theory is constructed so that all but one copy appear dark and gravitationally repulsive for any 
observer. We have shown that the cosmology of this multimetric theory naturally explains an 
accelerating universe. 

Motivation for our theory comes from the basic idea that Newtonian gravity allows positive 
and negative mass observers, but relativistic Einstein gravity does not (while it is well-known that 
it may contain negative mass sources). The main observation needed to transfer this idea into a 
general relativistic framework is that N different types of observers require for their definition N 
different metrics, as already argued in our article |25j . Corresponding to these different metrics 
one then also needs different copies of the standard model to distinguish the sources. Interestingly, 
the same article proves a no-go theorem that makes it impossible to construct bimetric gravity 
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theories in which a repulsive gravitational interaction between the different standard model sectors 
is contrasted by an attractive gravitational interaction of equal strength within each sector. Hence 
the first important result of this paper is the proof, by explicit construction, that this no-go theorem 
cannot be extended beyond N = 2, simply because canonical extensions of Einstein gravity with 
the stated properties exist for N > 3. 

This construction, as far as we are aware, provides the first relativistic gravity theory with local 
repulsive forces and suitably reacting observers which does not contain additional fields without 
clear interpretation, but merely extra dark copies of the well-understood standard model. Thus 
the dark universe is constituted by the same type of matter known from the visible universe, rather 
than that the visible universe is distinguished from the dark side by its physical properties. 

Since repulsive gravity is contained in our theory one could expect it to explain dark energy. 
To see whether this is indeed the case, we have analyzed a simple homogeneous and isotropic 
cosmological model. We have argued that the very early and the very late universe which are 
in all sectors dominated by radiation or matter, respectively, can be plausibly described by a 
single effective metric and energy momentum tensor. Essential ingredient of this argument was an 
extended Copernican principle of symmetry between all sectors with regard to initial conditions. 
We have derived the cosmological equations of motion and computed the general solutions for 
radiation and dust matter. Both solutions feature an open universe and the scale factor displays a 
big bounce. But most remarkably, this simple model also leads to accelerating expansion and the 
acceleration becomes naturally small at late times. This effect does not rely on exotic new matter; 
it is a direct consequence of the mutual repulsion between the different standard model copies in 
the theory. 

In this article we have discussed only the simple setting of a homogeneous and isotropic universe 
whose density and pressure functions effectively agree within each sector. Since this is plausible 
only for the early and late universe one should also consider the more general cosmological setting 
where these functions evolve independently. This requires a study not only of the simplified ef- 



fective equations (20), but of the more complicated full equations of motion (16) with the general 



cosmological ansatz (21) and (22). 

In further research, it would be interesting to apply the theory presented here to other im- 
portant astronomical observations on non-cosmological scales that involve non-homogeneous mass 
distributions. In this context one should discuss the consequences of the additional dark standard 
model copies for the physics of the solar system, for rotational curves of galaxies, gravitational 
lensing, or structure formation. Could it be that effects conventionally attributed to dark matter 
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simply follow from the dark standard model sectors? We know from the symmetry assumption 
between all sectors that all types of matter have the same physical properties, and thus should 
form the same structures, like stars and galaxies, as does visible matter. What we do not yet know 
can be learned by applying the comprehensive parametrized post-Newtonian formalism [28] to our 
theory. In our case, this formalism will involve two sets of PPN parameters: the first will describe 
the usual gravitational interaction within each sector, the second will specify the cross-interaction 
between different sectors. Once these parameters are available we will be in the position to link 
them to astronomical observations. 
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